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Abstract 


A mathematical model containing the essential features embod- 
ied in the noise suppression of lined ejectors is presented. Although 
some simplification of the physics is necessary to render the model 
mathematically tractable, the current model is the most versatile and 
technologically advanced at the current time. A system of linearized 
equations and the boundary conditions governing the sound field 
are derived starting from the equations of fluid dynamics. A non- 
reflecting boundary condition is developed. In view of the complex na- 
ture of the equations, a parametric study requires the use of numerical 
techniques and modern computers. A finite element algorithm that 
solves the differential equations coupled with the boundary condition 
is then introduced. The numerical method results in a matrix equa- 
tion with several hundred thousand degrees of freedom that is solved 
efficiently on a supercomputer. The model is validated by comparing 
results either with exact solutions or with approximate solutions from 
other works. In each case, excellent correlations are obtained. The 
usefulness of the model as an optimization tool and the importance 
of variable impedance liners as a mechanism for achieving broadband 
suppression within a lined ejector are demonstrated. 

Introduction 

The number of airplane passengers traveling in the long-range international market is ex- 
pected to quadruple within 25 years. International travel experts agree that an environmentally 
acceptable and economically viable supersonic transport could obtain nearly 30 percent of this 
lucrative market. This forecast has given development of a high-speed civil transport (HSCT) 
a new impetus within the United States, Europe, and the Orient. Such an airplane could be 
expected to fly as soon as the early part of the next century and to cruise at speeds between 
Mach 2 and Mach 3. The first country to put such a technology in place would be in a leader- 
ship position in long-range civil aviation for many years to come. NASA and several private 
companies have taken up the challenge of putting such a technology in place for the United 
States. The current focus of NASA is on ozone depletion and noise; the latter is the subject of 
this paper. 

To become environmentally acceptable, the HSCT must achieve noise levels comparable to 
those of other large commercial airplanes expected to fly during that period. Several studies 
(refs. 1-5) show that a reduction of 20 decibels in the radiated noise levels on takeoff will 
be required before the HSCT will meet stage-III noise standards. The mixer-ejector exhaust 
component is expected to provide one means of reducing the jet noise by mixing the high- 
velocity gases of the exhaust jet with much cooler ambient air; the result is a stream of air with 
a lower noise and energy level. The mixer-ejector is expected to provide only about 12 of the 
20-dB reduction requirement. (See refs. 6-9.) Thus, a second component, which is expected to 
be a lined ejector, will be required to bring the radiated noise to environmentally acceptable 
levels. 

Experimental tests, which would be required to develop an acceptable lined ejector technol- 
ogy, are expensive in terms of both research dollars and time. Accurate mathematical mod- 
els for predicting and minimizing the noise field within a lined ejector can be important in 
this development. These models could be used as inexpensive, nonlabor-intensive technology 
tools to optimize the lining for maximum noise suppression, to explore more innovative lining 
concepts, and to minimize expensive experimental tests. Mathematical models currently exist 
for predicting liner noise suppression properties in flow; these models are restricted mostly to 



irrotational subsonic flow and to sources that are tonal in nature such as fan noise. Several arti- 
cles representing only a portion of what is available have been cited in references 10-15. Present 
mathematical models are not applicable for the analysis of the lined ejector noise suppression 
because they do not properly account for the characteristics of the jet exhaust that will generate 
this complex noise field. 

The exhaust jet of the HSCT is expected to flow within a rectangular domain and will 
consist of an extremely hot supersonic core superimposed above and below with a cooler subsonic 
flow. The high-frequency broadband noise generated by the intense mixing of this composite 
flow field within the ejector escapes to the far field. This portion of the noise spectrum is 
targeted for reduction by the lined ejector. Aspects of liner noise reduction that are new to 
this environment include a rotational composite flow field and an extremely hot flowing gas. 
In addition, the jet noise source is distributed, broadbanded, and more directed than noise 
sources considered to date. Variable impedance liners, as considered in this work, are expected 
to be an essential component for attenuation of broadband jet noise. (See refs. 16 and 17.) An 
intense literature search produced no mathematical models for predicting the attenuation in this 
expected environment. Note that in the 1970’s, Abrahamson (refs. 18 and 19) developed the 
most promising mathematical model for analyzing rotational flows in acoustically treated ducts 
with variable impedance liners. However, this model is restricted to subsonic flow and circular 
axisymmetric geometries, tailored to fan duct noise sources, and developed for a computer system 
that is more than 25 years old. Furthermore, recent analyses (ref. 20) show that an additional 
boundary condition is required for accurate solutions to vortical flow problems. 

This paper describes a new mathematical model that simulates the essential characteristics 
of noise suppression with lined ejectors. A system of equations governing the sound field in a 
viscous, composite, vortical flow field is derived starting from the equations of fluid dynamics. 
Boundary conditions (including a new nonreflecting condition) consistent with this system are 
then introduced. A numerical method solves the differential equations with boundary conditions 
on a supercomputer. Several example problems are modeled and the results compared with their 
exact analytical solutions; comparisons with other approximate solution methods are made when 
exact solutions are not possible. Supersonic flow analysis is validated by comparing the model 
results with the results of a recently developed modal theory for plug flow. (See ref. 21.) Finally, 
the usefulness of the model as an optimization tool is tested. The model is also used to show 
the importance of variable impedance liners as a mechanism for achieving broadband noise 
suppression within an ejector. 

Problem Formulation 

A typical jet exhaust enshrouded by a lined ejector is illustrated in figure 1. The fluid 
within the mixing duct is considered to be viscous, composite, laminar, and compressible. The 
flow field within the ejector is considered a composite of subsonic and supersonic flow. To 
reduce computational requirements for the solution, only two-dimensional flow is considered; no 
variations in the flow occur in the direction normal to the X-Y plane. The inflow and outflow 
boundaries are at x — 0 and x — L, respectively. Sound sources occurring within the region 
0 < x < L will be modeled accordingly. The lower and upper boundaries at y — 0 and y — H , 
respectively, contain the wall lining. Noise suppression properties of the lining are specified by 
the normal acoustic impedances of the upper and lower walls, zjj and zq, respectively. These 
impedances are functions of the frequency content of the acoustic waves, position along the 
wall, and steady-state flow Mach number. In this paper, a numerical method is described for 
predicting the sound attenuation within the lined region 0 < x < L and for optimizing the wall 
lining. 
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Governing Equations 


The two-dimensional Navier-Stokes equations which govern the fluid flow within the domain 
described in figure 1 are 


d{U} d{H} d{G} 
dt dx dy 


= { 0 } 


(1) 


The vectors in equation (1) are defined elsewhere (e.g., ref. 22) and are not repeated here. The 
unknown variables in {U} are the fluid density p, axial component of velocity u, transverse 
component of velocity v, and total energy E . 


To derive the equations governing acoustic disturbances in the total flow field, the unknown 
fluid variables are resolved into the following steady flow field and steady-state acoustic field 
components: 
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where vectors {$}, {fl>o}, and {$} contain the total flow field, steady flow field, and steady- 
state acoustic field in the flow, respectively. Here, / is the frequency in hertz and i — — 1. The 

following observations are made concerning equations (2): the steady flow field is a function of 
both spatial coordinates and the normal assumptions for parallel and/or irrotational duct flow 
are not made; the axially propagating wave assumption is not assumed, which means that vortical 
acoustic disturbances are possible; and an ideal gas is assumed which eliminates temperature as 
a variable. 


After substituting equations (2) in equation (1), the following linearized system governing 
the acoustic disturbances in the flow field is obtained: 
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The ambient speed of sound cq is defined by 
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Equations (3) may be interpreted as the linearized steady-state counterpart of that derived in 
reference 20. 

The effects of nonlinearity, viscosity, and heat conduction are important in establishing the 
steady flow field and sources of sound, but they are not of primary importance in determining 
the steady-state acoustic disturbances discussed in this paper. Although some solutions of 
equation (1) for a general initial condition contain both the steady-state and the transient wave 
fields, viscous forces present in equation (1) will rapidly dissipate the transient wave field and 
leave only the steady-state acoustic field of equations (3) as the dominant noise field. Note that 
determination of the steady flow field {$o} is necessary for solution of the acoustic problem, but 
the assumptions and techniques used to obtain this solution are not part of this work. Instead, 
the steady flow field is assumed to be available in a form suitable for inclusion into this equation. 
Finally, before equation (3a) can be solved, a set of boundary conditions consistent with this 
system must be specified. 

The computational domain consists of two physical boundaries (the acoustically lined walls) 
and two artificial or free boundaries (the inflow and outflow boundary). The solution of the 
acoustic problem requires three boundary conditions along the subsonic portion of the inflow 
boundary and four along the supersonic portion. (See ref. 20.) Note that in references 18 
and 19, only two boundary conditions were imposed at the subsonic inflow boundary. In general, 
two boundary conditions are insufficient for acoustic calculations in vortical flows. Along the 
outflow boundary, no boundary condition is needed for the supersonic portion; however, a single 
boundary condition is required along the subsonic portion. 

As a result, along the inflow boundary acoustic pressure and tangential component of acoustic 
velocity are specified and the disturbance is assumed homentropic as follows: 


p(0,y) - p s {y) 

(5a) 

v(0, y) - v s (y) 

(5b) 

p(°>y) = c 2 0 (0,y)p(0,y) 

(5c) 


In equations (5a) and (5b), p s (y) and v s (y) are the source acoustic pressure and tangential 
component of acoustic velocity, respectively, at the inflow boundary; equation (5c) is the 
homentropic inflow condition. In addition to equations (5), two conditions are required: one 
along any supersonic portion of the inflow boundary and one along the subsonic part of the 
outflow boundary. These additional conditions are presented after the discussion of the numerical 
method. 

The acoustic boundary condition at the treated walls must be determined by the nature of 
the interaction of the acoustic waves and the mechanical structure forming the wall. Low-speed 
subsonic flow can possibly result in the ratio of acoustic pressure to the normal component of 
acoustic velocity at the surface assuming a specified value, which is referred to as the normal 
acoustic impedance. Whether this can be achieved for a composite flow containing an extremely 
hot flowing gas is currently a research issue. The current model assumes that a locally reacting 
surface impedance can be defined for the lining material even in this new environment. For no 
slip at the liner and flow interface, the condition can be stated as 


p(x,0) 

S(z,0) " ( * 

(6a) 
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nx,H) = ZHt - I) 

(6b) 
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Occasionally, a uniform or plug flow model will be convenient to permit fluid slip at the wall. 
The proper form of the wall boundary condition (ref. 23) for this situation is 


v(x, 0) 


p(x, 0) uo(£,0) d p(x, 0) 

zq(x) iu> dx\_ zq(x) 


(7a) 


v(x, H ) 


p(x,H ) uq(x,H ) d p(x,H) 

Zff(x) iio dx\_ Zff(x) 


(7b) 


where the component of the steady velocity uq in the Y direction is assumed zero at the wall. 
As expected, equations (7) reduce to equations (6) when there is no slip at the wall. 

Equations (3) describe linear acoustic disturbances in viscous, rotational, and compressible 
fluid flow. Note that in the absence of any disturbance, equation (3a) is satisfied identically. 
No exact analytical solutions of these equations exist for an arbitrary steady flow. The 
method chosen to solve the acoustic equations (3) coupled with boundary conditions defined 
by equations (4)-(7) is a Galerkin finite element method, which is discussed in the following 
section. 


Galerkin Finite Element Analysis 

Details on the finite element theory may be found in several texts. (See, for example, refs. 24 
and 25.) The finite element method, when applied to the current acoustic problem, represents 
the continuous sound field as an assemblage of rectangular elements interconnected at nodes 
as illustrated in figure 2(a): N nodes in the axial direction and M nodes in the transverse 
direction. Note that the acoustic field is zoned at imaginary lines although no physical separation 
is envisioned. A typical rectangular element e is shown in figure 2(b). Each element consists 
of four nodes labeled 1, 2, 3, and 4, respectively, and is considered to have width a and height 
6. Values of a and 6 for each element should be chosen so that the acoustic disturbances and 
steady flow field within each element are simulated as closely as possible. The objective is to 
obtain the unknown acoustic variables at the nodes of each of the (M — 1), ( N — 1) elements. 

Galerkin’s finite element method is used to minimize the error vector and reduces the problem 
to a finite set of algebraic equations which are solved using matrix methods. The error vector is 
defined as _ _ 

{M*,S)} = + [«0]fi + [C 0 ] {4} (8) 

Within each element, the acoustic pressure, two components of acoustic velocity, and the density 
will be represented as the following linear functions: 
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in which pj, uj, vj, and pj are the values of acoustic pressure, axial component of acoustic 
velocity, transverse component of acoustic velocity, and acoustic density at node J , respectively. 
The steady flow field is represented in a similar manner as 

4 4 

Po(x,y) - ^2 Nj(x,y)Pj u 0 (x,y)-^2Nj(x,y)Uj (9e) 

J = 1 J = 1 

4 4 

u 0 (x,y) - ^ Nj(x,y)Vj p 0 (x, y) = ^ Nj(x, y)Rj (9f) 

J = 1 J = 1 


The correct solution of equations (3) for the acoustic system is obtained when the error 
vector {E r (x, y)} is identically zero at each point of the domain. An approximate solution is 
achieved by requiring that the error vector be orthogonal to each basis function N j(x, y) that is 
assumed to represent a complete set of functions. The contribution of a typical element to the 
minimization of the error function is 

[ [ {E r (x,y)}Nj(x,y)dydx = [A e ]{<S> e } 

Jo Jo 

in which [ A e ] is a 16 X 16 complex matrix (i.e. , the stiffness matrix) and {<I> e } is a 16 X 1 
column vector containing the unknown acoustic variables at the four nodes of the element. The 
coefficients in the local stiffness matrix [A e ] were computed in closed form and are not included 
in this report. 

Assembly of the global stiffness matrix for the computational domain is a basic procedure in 
the finite element method. Appropriate shifting of rows and columns is all that is required to 
add the local element [A e ] directly into the global matrix [A]. Assembling the elements for the 
entire domain results in a matrix equation of the form 

= (io) 

in which [A] is a 4MN X 4MN complex matrix, and {$} is a 4MN X 1 column vector containing 
the nodal values of the unknown acoustic variables. In equation (10), the 4MN X 1 column 
vector {F} is identically zero and the global stiffness matrix [A] is singular. Therefore, boundary 
conditions must be applied to this system of equations before a solution can be obtained. 

As mentioned in the previous section, an additional condition must be specified before 
the solution can be calculated. Determining the inflow or outflow boundary condition is 
difficult. The difficulty arises because the sound field within the ejector propagates through 
the multidimensional steady flow field and radiates outward to the surrounding environment. 
The additional boundary condition must model this characteristic. Several methods have been 
proposed for determining this boundary condition in subsonic flows. Abrahamson (refs. 18 
and 19) expressed this boundary condition as a local linear relationship between the acoustic 
pressure and the normal component of acoustic velocity. The method of determining the constant 
in this linear relationship was not specified. Horowitz, Sigman, and Zinn (ref. 26) expressed this 
boundary condition in terms of an exit impedance. Because the exit impedance was not known 
a priori, the solution within the finite computational domain was combined with an integral 
formulation for the far field. The exit impedance was then compared with that obtained from 
the integral formulation at the boundary and was iteratively corrected until the results converged. 
A similar method was applied by Baumeister. (See ref. 27.) 
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In this report, the additional condition is expressed as a linear relationship between the 
acoustic pressure at a boundary node and the normal component of acoustic velocity along the 
boundary nodes: 



{PB} = [Z]{UB} 


(11a) 

[Z] = 

’ -Zn z i2 
Z 21 Z 22 

Z 1M " 
Z 2M 

(lib) 


-Z M \ %M2 ■■■ 

Z MM - 



Here {PB} and {U B} are 1 X M column vectors containing the acoustic pressure and the 
normal component of acoustic velocity, respectively, along the boundary nodes. Determining the 
complex coefficients in the MxM matrix [Z] is difficult for the general case. The elements in this 
matrix should result in the acoustic held having the correct radiation condition. The coefficients 
in [Z] that ensure no reflections are derived in the appendix. This derivation assumes that the 
duct terminates in an environment that is homogeneous in the axial direction. Alternatively, [Z] 
could be obtained by an extension of the iterative method proposed in reference 26. 

Numerical Implementation of Boundary Conditions 

Along the boundaries of the computational domain, three types of boundary conditions may 
exist: the noise source pressure or tangential velocity, the homentropic acoustic wave condition 
of equation (5c), or the additional condition of equations (11). All three boundary conditions 
can be inserted into the assembled global matrix equation (10) by imposing relations on the 
nodal values of acoustic variables along the boundary. Details for imposing the boundary 
conditions are described in references 18, 24, and 25 and are not repeated here. The noise 
source boundary condition consists simply of specifying all nodal values of acoustic pressure 
or velocity. The homentropic acoustic wave condition involves specifying a linear relationship 
between acoustic pressure and acoustic density at the inflow boundary nodes. The additional 
condition in equations (11) involves specifying a linear relationship between acoustic pressure at 
a node and the normal component of acoustic velocity at the boundary nodes. 

Solution to the Matrix Equation 

The global matrix [A], generated by Galerkin’s method after the application of boundary 
conditions, is an unsymmetric, positive indefinite, complex matrix. Fortunately, because of the 
discretization method used, the matrix will be square and block tridiagonal of order 4MN ; the 
structure of matrix [A] is shown in figure 3(a). This global matrix contains a number of major 
blocks aj, bj, and cj that are square and block tridiagonal and have the structure shown in 
figure 3(b). Each minor block Aj, Bj, and D j is a square matrix the order of which is four. 
Of practical significance, the structure minimizes computer storage requirements and maximizes 
computational efficiency. All computation, storage, and boundary condition manipulations are 
performed on the block triangular portion only of the matrix [A]. (See fig. 3(a).) Special matrix 
techniques exist for a solution of this structure. Gaussian elimination with partial pivoting 
and equivalent row infinity norm scaling is used to reduce the rectangular system to an upper 
triangular form. Back-substitution is then employed to obtain the solution for the acoustic 
variables. All computations were performed on the supercomputer. 

Evaluation of Liner Performance 

A primary objective of the numerical method is to evaluate the performance of the wall 
lining. The following uniform flow expression (ref. 28) is used for the acoustic intensity at a 
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point (x, y) in a uniform steady flow: 

I(x, y) - +p*u) + ^y~Poc o (uu* + vv * + (12) 

4 4 V p 0 c 0 J 

The steady flow Mach number Mq is defined as uq/cq and the superscript asterisk denotes the 
complex conjugate. The attenuation of the lining in decibels is then obtained from the equations 


AdB = 10k.g lo ^) ) 

(13a) 



W{x) - / I(x, y)dy 

(13b) 


Jo 


To obtain the axial acoustic energy flux W(x), the integration in equation (13b) was performed 
in closed form. Equations (13) are used to evaluate the effect of the wall lining in this report; 
results in the next section are purposely restricted to uniform flow. 

Results and Discussion 

To validate the numerical model, finite element method results are compared with exact 
analytical results for a rigid wall duct and with an approximate modal theory recently developed 
for uniform supersonic flow between soft walls. (See ref. 21.) These solutions are compared 
for a range of frequencies and Mach numbers. Such comparisons provide confidence in the 
boundary condition formulation, convergence characteristics, computational cost, and overall 
accuracy of the numerical model. The model was also used to perform optimization studies and 
to demonstrate that the variable impedance liner is a candidate for enhancing the broadband 
performance of lined ejectors. 


Results for a Rigid Wall Duct 

For a rigid wall duct (zq(x) — Zff(x) — ► oo) with uniform flow in the positive X direction, 
the acoustic pressure field consisting only of a simple progressive single mode is 


p(x, y) - A m cos(^p)e* A 


(14) 


where K m satisfies the condition 


X m — 


-kM + x jk 2 M 2 + {l- M 2 ) k 2 -( Vff) 
M 2 - 1 


and 


M — — k = — 


(m — 0,1,2,..., oo) 


co c 0 

Equation (14) may be regarded as a special solution of equations (3) with 


Ps(y) - T m cos(^p) P-clp v s (y) 


ip'siv) 


Po(ui + K m u 0 ) 

The matrix \Z], which ensures the existence of a simple progressive traveling wave, is 


W = 


Po(lv + K m u q) 
X m 


[I] 
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where [/] is the identity matrix and only a single mode to is assumed to propagate. The above is 
a standard result that is based on a constant steady flow field and is computed by separation of 
variables. The solution is useful for comparison with numerical calculations. Numerical results 
for a rigid wall duct were computed with ambient values of po, Po, and 2~b, with L — H — 0.5 m 
and A m — 1; and with an evenly spaced 51 x51 mesh. Results were computed for a subsonic and 
supersonic Mach number (Mg = 0.5 and Mg = 1.5), three frequencies (/ = 500 Hz, / = 5000 Hz, 
and / = 10 000 Hz), and a planar (to = 0) and nonplanar (to = 1) source. Although several 
transverse and axial locations were checked at random in the duct, results are shown only for 
the transverse location y — 0.75 H to limit the number of graphs. 

Graphs of the real and imaginary solutions for the planar (to = 0) acoustic pressure computed 
from this analysis are compared with the exact steady state solution for the subsonic and 
supersonic Mach number in figures 4 and 5, respectively. Computations are for a frequency 
of 500 Hz and excellent correlation was obtained between the exact and numerical solution. No 
reflections as a result of the chosen boundary condition are evident in the figures. However, note 
that the wavelength of the acoustic wave in the supersonic flow has been stretched compared 
with that of the subsonic flow. Thus, the number of computational points required for accurate 
resolution may be considerably less for sound waves in a supersonic flow. 

Results are shown in figures 6 and 7 when the frequency of the planar wave is increased an 
order of magnitude (/ = 5000 Hz). Again, excellent correlation was obtained for the subsonic 
and supersonic Mach number. Note that for the subsonic Mach number, approximately 10 axial 
points per wavelength (i.e. , N — 51) are used in the finite element computer code. The 10 axial 
points per wavelength are considered the least number of axial nodes for which an accurate 
numerical solution can be expected in the subsonic flow. Figure 8 shows the results for the 
supersonic Mach number at the highest frequency (/ = 10 000 Hz). Good correlation is again 
obtained. Numerical results were also obtained for a nonplanar (to = 1) wave source at all three 
frequencies. Figure 9 shows a sample calculation for the nonplanar wave at the highest frequency 
(/ = 10 000 Hz) and the supersonic Mach number. Good correlation between the finite element 
and exact solutions is obtained; better correlation than that of figure 9 was obtained at the two 
lower frequencies, although results are not shown for the sake of brevity. 


Results With Wall Lining 

Lining materials used in the ejector of the HSCT should have the following characteristics: 
have minimum thickness and weight, resist high pressures, withstand extreme temperatures, 
and absorb broadband sound. Results in the previous section confirmed that the mathematical 
model can be used to calculate in detail the sound field for planar and nonplanar acoustic waves. 
In this section, the application of the model for prediction of the sound attenuation of a wall 
lining is discussed. Although no exact solutions are possible for vortical flow, an approximate 
modal solution is possible in uniform flow. This modal solution has been explored in detail by 
several other researchers for subsonic flows. Recently, the modal solution (ref. 21) was expanded 
to include uniform supersonic flow. This modal solution is used as a basis to test predictions of 
the model discussed in this report. Because data will soon be available from ejector tests to be 
performed in the Langley Jet Noise Facility, the geometry of that facility (i.e., L — 0.3302 m 
and H — 0.127 m) was used to develop the results that follow. Further, a 51 X 51 evenly spaced 
grid was used in all calculations. 

A comparison of the finite element model attenuation predictions with the modal theory was 
made for a planewave source (i.e., p s (y) — 1) at 1500 Hz with zjj — zq — 2 — 0.58*. The modal 
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theory represented the source pressure with a four-mode expansion. Comparisons were made by 
using the following [Z] matrix: 

zi 0 ... 0 

0 z 2 • • • 0 

0 0 . . . z M . 

in which 

P(x, Vi) 

u(x, t/i) 

Values of p(x, y z ) and u(x,yj) were computed from the four-mode expansion. The matrix \Z] 
above has the advantage that only four modes in the lined section have to be computed. The 
disadvantage of this choice of [Z] is that it requires prior knowledge of the solution. Note that 
although the full matrix [Z] as derived in the appendix is preferred, the 51 soft wall modes 
required for the computation were not readily available from the modal theory. 

Results are plotted as a function of Mach number in figure 10. As shown, the attenuation 
is reduced with increasing Mach number even in the supersonic flow. Again, note that liner 
attenuation in supersonic flow as considered here represents a new development. Further, the 
finite element and modal theory show no discrepancy in the subsonic flow. Both predictions 
correlate closely in the supersonic flow region, although the small oscillation that occurs in the 
finite element results is not characteristic of the modal theory. Similar characteristics were 
observed for most of the other frequency values. 

An Optimization Study on Broadband Suppression 

The HSCT requires that special attention be given to optimization of liner properties. Unlike 
fan noise sources that are tonal in nature, jet sources are distributed and broadband. The 
impedance of a lined ejector must be determined so that maximum broadband suppression is 
achieved. The present model is applicable to such an extensive and thorough parametric study. 
To demonstrate the mathematical model for this purpose, an analysis to optimize a uniform 
lining impedance was completed. An optimal uniform liner has been obtained by mapping 
contours in the complex impedance plane. A segmented lining was then developed which had 
greater broadband suppression than that of an optimal uniform lining. 

Figure 11(a) shows a variable impedance liner designed to enhance the broadband perfor- 
mance of the ejector. For simplicity, a two-segmented liner with each segment 0.5 L units long 
has been chosen. Channel dimensions (L and H ) were equivalent to those of the ejector located 
in the Langley Jet Noise Laboratory. The impedances of the first and second segment of the 
segmented liner were not optimized but were taken from the finite element model results of the 
individual uniform liner optimization studies in a supersonic flow. 

Figure 11(b) shows the result of the design at M o = 1.5 and frequencies of 0-10 000 Hz. A po c O 
entrance impedance and a planar source pressure have been assumed. The curve data for the two 
uniform liners were computed by giving each segment in figure 11(a) the same wall impedance. 
The uniform liner was then optimized for maximum suppression at a fixed frequency defined 
as the tuning frequency. Both optimized uniform liners were observed to perform well at their 
tuning frequency but their performances fall off rapidly on either side of their tuning frequency. 
This characteristic is clearly undesirable for the broadband jet noise source which would occur 
in an ejector. Attenuation predictions for the segmented liner with the first segment tuned at 
2500 Hz and the second segment tuned at 5000 Hz are also shown in the figure. The segmented 
liner suppresses more high-frequency sound than the uniform liner tuned at 2500 Hz and more 
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low-frequency sound than the uniform liner tuned at 5000 Hz. Therefore, the segmented liner 
is more desirable for a broadband jet noise source. If the impedance of the first and second 
segment were actually optimized, the attenuation spectrum for the segmented liner could be 
much broader. 

Conclusions 

The objective of this report was to introduce the equations and develop the numerical method 
for predicting and optimizing the noise suppression of lined ejectors. A complete validation of the 
full capability of the model in a single report is not possible; the results were purposely restricted 
to irrotational flows because exact and approximate solutions are available in the literature for 
comparison and an expression for the liner performance is known. Validation of the model in a 
vortical flow field has not yet been accomplished and is beyond the scope of this report. Results 
from the model to date are satisfactory and show expected convergence characteristics over a 
wide range of frequencies, Mach numbers, and source conditions. 

The following specific conclusions are based upon these results: 

1. A mathematical model believed to contain the essential features required for predicting sound 
suppression by lined ejectors is developed. 

2. The model is compatible with a composite flow field (vortical and nonvortical), variable 
impedance liners, arbitrary noise sources, and the current generation of supercomputers. 

3. Although a simplification of the physics has been necessary to make the problem mathemat- 
ically tractable, the current model is the most versatile and technologically advanced at the 
present time. 

4. Excellent correlation with the exact solution for planar and nonplanar sources in subsonic and 
supersonic flow in a rigid wall duct is obtained using approximately 10 points per wavelength. 

5. Results with wall lining show excellent correlation with a recently developed modal theory 
(ref. 21) for uniform supersonic flow. 

6. The use of the model as an optimization tool is successfully demonstrated. 

7. Optimization studies with the model confirm the advantage of variable impedance liners as 
a possible mechanism for achieving broadband noise suppression in lined ejectors. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
February 17, 1994 
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Appendix 

Nonreflecting Condition for Inflow and Outflow Boundaries 

In this appendix, elements for the coefficient matrix [Z] in equations (11) of the body of this 
report are derived. Boundary conditions assume that the transverse component of the mean 
velocity ro is identically zero and p q, To, and p q are constant. The mean flow uq is considered a 
function of the transverse coordinate only and the disturbance is assumed homentropic. Finally, 
the boundary is assumed to terminate in an atmosphere that is homogeneous in the axial 
direction. Under these assumptions, the solutions to equations (3) in the form of outgoing 
waves are of the following complex exponential form: 

P(x,y)= '£, A ™^T^- eiKmX (Al) 

m= 1 C 0 

M 

u(x, y) - ^2 B m Um(y)e lKrnX (A2) 

m= 1 
M 

v(x,y)=Y / D mV m (y)e iKmX (A3) 

m= 1 
M 

P(x, y) - E A mPm(y)e lKmX (A4) 

m= 1 

Here, K m is the axial propagation constant and U m {y), V m {y), and P m (y) are the eigenfunctions 
associated with the axial component of acoustic velocity, transverse component of acoustic 
velocity, and acoustic pressure, respectively. Note that each series has been truncated at a 
finite number M, where M is the number of nodes along the boundary. To ensure no reflections, 
the sums in equations (A1)-(A4) are taken only over modes whose axial propagation constant 
possesses positive imaginary parts. 

A nonreflecting condition relating the acoustic pressure and axial velocity at a boundary 
node is now developed from these series expansions. Substitution of equations (A1)-(A4) into 
the second component of equations (3) and elimination of v give 


M 

El i[u> + K t 

m = 1 


iuo(y)\B m Um T 


iPm A m,Pm 


P 0 


u o (y) A mPUy) \ _ Q 
ipo[u) + K mM o(y)]J 


(A5) 


Multiplying equation (A5) by each P n {y) and integrating across the boundary give 


/0 


H f T , w ( m d tt , iPm A mPm{y) U Q {y) A mPm{y) l D MJ n f \ a\ 

}_ J ) + h m u 0 (y)\B m U m H : ^ } Pn{y ) ay — 0 (A6) 


m= 1 


P 0 


ip 0 [w + K rn U 0 (y)] 


Equation (A6) consists of M equations in 2 M unknowns that may be expressed in matrix format 
as 


{A} = [Z b ]{B) 


{A}={ 


A i 

A 2 


{B}= { 


Bi 

b 2 


{ a M 


Bm 


(A7) 
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- a n 

«12 • 

• • a lM " 

-1 

' &n 

&12 

■ ■ b 1M ' 

[Zb\ = ~ 

«21 

«22 

a 2M 


hi 

^22 

• • b 2M 


-“Ml 

“M2 • 

•• “MM- 


- b M 1 

^M2 ■ 

•• ^MM - 


(A8) 


f H 

l mn — / l[u + K m U 0 (y)]Um{y)Pn{y) dy 
JO 


rH 


iK m , s u'niv) 

P m (y )-- ° ’ 


I o l PO ' ipo[u) + K m u 0 (y)Y 
The vectors {PB} and {UB} in equations (11) are 

{PB} = [P][E]{A} 
{UB}=[U][E]{B} 


P'm(y) \ Pn(y)dy 


[E] 


giK-yX Q 

q e iK 2 x 


[P\ = 


[U] 


. o 

0 

e iK M x J 

" -Pi(yi) 

P2(yi) 

• • • Pm{vi) ' 

■Pi(y 2 ) 

P2iV2) 

■ ■ ■ Pm{V2) 

-Piiyu) 

P2{um) 

••• Pm{1)m) - 

' Ui(yi) 

U2iVl) 

■ ■ ■ U M (yi) " 

Uiim) 

^2(1/2) 

••• P Aiiy 2 ) 

U\{yM) 

U2{yM) 

• •• Ua {( yu). 


(A9) 

(A10) 

(All) 

(A12) 

(A13) 


(A14) 


(A15) 


in which the point y — yj corresponds to a boundary node point and x corresponds to the axial 
location of the boundary. Substitution of equation (A7) into equation (All) gives 


{PB} = [P][E][Z b ]{B} 

Solution of equation (A12) for {B} and substitution into equation (A16) give 


(A16) 


{PB} = [. Z]{UB } (A17) 

[Z] = [PHEpiM-’p]- 1 (A18) 

Equations (A17) and (A18) are the nonreflecting boundary conditions referred to in the body 
of this report. This general condition is not dependent upon the sound source. The matrix [Z] 
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can be determined provided that each of the M eigenfunctions P m (y), U m (y), and axial wave 
numbers K m is known. 

The equations governing these eigenfunctions are obtained by substituting the axially 
propagating wave solution into equations (3) and (7) to obtain 


+ h m uo(y)]Um+ +u 0 (y)V m -0 

P 0 

(A19) 

T A m t<o(y)].P m + p qc q(? A m U m + Vffi) — 0 

(A20) 

p' 

v - m 

m ipo[ui + K m u 0 (y)] 

(A21) 

v m {o) — ~ l [i + Uo(0)A -k(o) 

(A22) 

V m (H) = k f 1 + MH)Km \ Pm(H) 
Z H l LO J 

(A23) 


The first order system governed by equations (A19)-(A23) may be integrated numerically to 
obtain the eigenfunctions and axial wave number required to construct \Z\. 

As an alternative to solving the first order system of equations (A19)-(A23), the equations 
may be combined into an equivalent second order system for the acoustic pressure eigenfunction 

P ™ iy) + p ' miy) + {fc2 “ K ™\ 1 ~ M ° 2(y) l + 2kK ^ M o(y)}Pm(y) = o (A24) 

pL( o) = 1 + 2 p -(°) ( A25 ) 

PUH) = 1 + ^ Pm(H) (A26) 

where the mean flow Mach number M 0 (y) and free space wave number k is defined by the 
equations 

k = - M 0 (y) = (A27) 

co c 0 

This second order system may be integrated numerically to obtain the pressure eigenfunction 
Pm(y) and axial decay rate K m . Equations (A20) and (A21) are then solved to obtain the 
eigenfunctions U m {y) and V m (y). 

The derivation of [Z\ is of interest for the special case of a uniform mean flow for which uq 
is a constant along the boundary. Because the mean flow gradient term u'q vanishes, 

D ( ^ r-M- + K m u 0 )] 

Pm(y ) = p Um{y) 

bmn — ~ a mn 
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Thus, [Zb] — [/] where [I] is the identity matrix and [P] and [U] are related through 
matrix [Z u ] 

[P] = [U][Z U ] 


[Zu] = 


Po( u +Ki u o) 

Kl 

0 


0 

Pq(lo+K 2 uq) 

k 2 


0 0 

For the special case of uniform flow, equation (A18) gives 


0 

0 


P 0 (lo+K m u o) 

k m 


[Z] = [P][C]-‘ 


Elimination of [U] from equation (A30) gives 

[z] = [PKZ.KP ]- 1 

Equation (A31) is the matrix [Z] for no reflections in a uniform flow. 


diagonal 

(A28) 

(A29) 

(A30) 

(A31) 
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(b) Each major block in matrix [A]. 
Figure 3. Block tridiagonal structure. 
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Imaginary part of acoustic pressure, Pa Real part of acoustic pressure, Pa 



(a) Real solution. 



(b) Imaginary solution. 

Figure 4. Comparison of planar acoustic pressure at 500 Hz and Mq — 0.5. 
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Exact solution 
Finite element solution 



Exact solution 
Finite element solution 



Axial coordinate x, m 

(b) Imaginary solution. 

Figure 8. Comparison of planar acoustic pressure at 10 000 Hz and Mq — 1.5. 
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0 .5 1.0 1.5 2.0 2.5 3.0 

Mach number, M 


Figure 10. Comparison of liner attenuation of a planar noise source at 1500 Hz with outgoing wave termination 
impedance. 
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Segmented liner ejector 



Goal: 12 dB (enhanced mixing) 

+8 dB (liner absorption) 

20 dB (suppression) 

(a) Variable impedance model cross section. 



(b) Attenuation with two-segment unoptimized liner; Mq — 1.5. 
Figure 11. Segmented liner for improved broadband performance. 
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